Evolution from BCS to BEC superfluidity in p-wave Fermi gases 
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We consider the evolution of superfluid properties of a three dimensional p-wave Fermi gas from 
weak (BCS) to strong (BEC) couphng as a function of scattering volume. We analyse the order 
parameter, quasi-particle excitation spectrum, chemical potential, average Cooper pair size and the 
momentum distribution in the ground state (T = 0). We also discuss the critical temperature Tc, 
chemical potential and number of unbound, scattering and bound fermions in the normal state 
(r = Tc). Lastly, we derive the time-dependent Ginzburg-Landau equation for T ~ Tc and extract 
the Ginzburg-Landau coherence length. 



PACS numbers: 03.75.Ss, 03.75.Hh, 05.30.Fk 



Arguably the next frontier of research in ultracold 
Fermi systems is the search for superfluidity in higher 
angular momentum states. Substantial experimental 
progress has been made recently Q, S S 13 in con- 
nection to p-wave cold Fermi gases, making them ideal 
candidates for the observation of novel triplet superfluid 
phases. These phases may be present not only in atomic, 
but also in nuclear (pairing in nuclei) , astrophysics (neu- 
tron stars), and condensed matter (organic superconduc- 
tors) systems. 

The tuning of p-wave interactions in ultracold Fermi 
gases was initially explored via p-wave Feshbach reso- 
nances in trap geometries for '^'^K P, Q and ^Li 0, 0|. 
Finding and sweeping through these resonances is diffi- 
cult since they are much narrower than the s-wave {£ ~ 0) 
case, because atoms interacting via higher angular mo- 
mentum channels (£ ^ 0) have to tunnel through a cen- 
trifugal barrier to couple to the bound state Fur- 
thermore, while losses due to two-body dipolar [3, |^ or 
three-body P, |3| processes challenged earlier p-wave ex- 
periments, these losses were still present but were less 
dramatic in the very recent optical lattice experiment in- 
volving and p-wave Feshbach resonances 0. 

For a dilute Fermi gas, the magnetic dipole- 

dipole interactions between valence electrons split p-wave 
{i = 1) Feshbach resonances that belong to different me 
states 0. Therefore, the ground state is highly depen- 
dent on the detuning and separation of these resonances, 
and possible p-wave superfluid phases can be studied 
from the Bardeen-Cooper-Schrieffer (BCS) to the Bose- 
Einstein condensation (BEC) regime. For instance, it 
has been proposed j3, for sufficiently large splittings 
that pairing occurs only in me = and does not occur 
in me = ±1 state, while for small splittings, pairing oc- 
curs via a linear combination of the me = and me = ±1 
states. Thus, these resonances may be tuned and studied 
independently if the splitting is large enough in compar- 
ison to the experimental resolution. 

The BCS to BEC evolution in p-wave systems was 
recently discussed at T = for a two-hyperfine state 
(THS) 121] in three dimensions (3D), and for a single- 



hyperfine state (SHS) [H [13 

in two dimensions, us- 
ing fermion-only models. Furthermore, fermion-boson 
models were proposed to describe p-wave superfluidity 
at zero 0, 1^ and flnite temperature |l2l | in three dimen- 
sions. Unlike the previous models, we present a zero and 
flnite temperature analysis of SHS p-wave Fermi gas in 
3D within a fermion-only description, where molecules 
naturally appear as bound states of two-fermions. 

The Hamiltonian for a dilute SHS p-wave Fermi gas in 
3D is given by 
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where the pseudo-spin t labels the hyperfine state rep- 
resented by the creation operator aj^. , and b]^^ = 

Here, f(k) = e(k) — ^, where e(k) = 
/ (2M) is the energy of the fermions and /i is the chem- 
ical potential. The attractive interaction can be writ- 
ten in a separable form as T^(k, k') = — 47r(7r*(k)r(k') 
where g > 0. The function r(k) = rk(fc)ra(k) is a 
symmetry factor where rk(fc) = (fcA:o)/(fc^ -I- fcp), and 
ra(k) = li.o(k) is the angular dependence. In addition, 
ko ^ Ro^y where i?o the interaction range in real space, 
sets the momentum scale. Furthermore, the diluteness 
condition {tiRq <C 1) requires {ko/k-pY ^ 1: where n is 
the density of atoms and k-p is the Fermi momentum. 

In the imaginary-time functional integral formalism 
(?i = fee = 1 and /3 = l/T), the partition function can 
be written as Z = J D[a\ a\e^^ with an action given by 

S = i^dr [EkaM(^)(^r)ak,T(T) + ^(t)] . We first in- 
troduce the Nambu spinor ?A^(p) — (Op ^ a_p^|), and de- 
fine p — (k. We) to denote both momentum and fermionic 
Matsubara frequency W£ = {2t+l)n / (3. Furthermore, we 
use the Hubbard-Stratonovich transformation to decou- 
ple fermionic and bosonic degrees of freedoms. Then, 
we integrate over the fermionic part, and rewrite the 
bosonic field as a combination of r-independent Ao and 
T-dependent h.{q). Here, q = (q, Wf) with bosonic Mat- 
subara frequency vi = 2i?7r//3. 
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Performing an expansion of S to quadratic order in 
A(g), we obtain 



where the vector A.'^iq) is such that A.^q) = 
[At(g), A(— g)], and the matrix F^-^{q) is the in- 
verse fluctuation propagator. Here, Sq is the sad- 
dle point action given by 5*0 = P\Ao\^ / {8ng) + 
Ep [/3C(k)/2-Trln(/3Go"V2)] , where the inverse 
Nambu propagator is Gq^ = iwiao — ^(k)(T3 -I- 
A5r(k)(T_ -I- r*(k)Aocr+. The fluctuation term in the 
action leads to a correction to the thermodynamic po- 
tential, which can be written as rJgauss = f^o + ^^fluct with 
no = So/ 13 and ^ifluct = ^ ^ Eg lndet[F-i(g)/(2/3)]. 

The saddle point condition (55o/(5Aq = leads to an 
equation for the order parameter 



— = E 



__tanh^— , 



(3) 



where -^(k) = (C^(k) + |A(k)p)3 is the quasi-particle 
energy, and A(k) — Aor(k) is the order parameter. The 
scattering amplitude within a T-matrix formulation is 
f{k) — fc^/(— l/flp -I- Tpfc^ — ik^) for the p-wave channel, 
where Op is the scattering volume, and r-p has dimensions 
of inverse length. Using /(fc), we can elliminate g in favor 



of Op via the relation 
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where V is the volume. 

The order parameter equation has to be solved self- 
consistently with the number equation N — —dft/dfi 
which leads to two contributions to the number equa- 
tion N — Nq + iVfluct- -^0 = —d^o/d^ is the saddle 
point number equation given by 



^ no(k); no(k) ^ - /ff x tanh 



/3£;(k) 



2E{k) 



(5) 



where rio(k) is the momentum distribution. 
Similarly, A'fiuct — — 5r2fiuct/3/i is the fluc- 
tuation number equation given by iVfluct = 
-r ' Eg{5[det ¥-\q)]/dy}/ det ¥-\q). 

For T « 0, A^fluct is small (cx T^) compared to A^o ^ 
for any interaction strength leading to A^ « A'o. In 
Fig. we plot Ar — Ao/ep and /ir = m/cf at T = 
as a function of l/(fcpap), where ep = kp/{2M) is the 
Fermi energy. Here, we choose ko ~ 200A:f. Notice 
that the BCS to BEC evolution range in l/(fcpap) is 
~ ko/kp. The weak coupling /i = ep changes contin- 
uously to the strong coupling ^ = — l/(Mfcoap) when 
fcgOp ^1. In strong coupling, Up has to be larger 
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FIG. 1: Plots of reduced a) order parameter amplitude = 
Ao/ep and chemical potential /ir — ^i/tF, and b) average 
Cooper pair size fcp^pair at T = and GL coherence length 
fcp^QL at r = Tc in a logarithmic scale versus l/(fc|ap). 



than Op > 2/fcg for the order parameter equation to 
have a solution with /i < 0, which reflects the Pauli 
exclusion principle. In addition, the weak coupling 
Ao = 24(A:o/A:p)ep exp[-8/3 + ^A:o/(4fcp) - 7r/(2fc||ap|)] 
evolves continuously to a constant Aq = 8ep[eo/(9ep)]^/^ 
in strong coupling, where eo = /cq/(2M). The evolution 
of Aq and fi are qualitatively similar to recent T = 
results for THS fermion ^] and SHS fermion-boson Q 
models. Due to the angular dependence of A(k), the 
quasi-particle spectrum £'(k) is gapless (mini?(k) = 0) 
for /i > 0, and fully gapped (mini?(k) — for /i < 0. 
Furthermore, both Aq and /i are non-analytic exactly 
when fj, crosses the bottom of the fermion energy band 
/i = at l/{kpap) « 0.5. The non-analyticity does not 
occur in the first derivative of Aq or fi as it is the case 
in 2D IQJ, but occurs in the second and higher deriva- 
tives. Therefore, the evolution from BCS to BEC is not 
smooth, and a top olog ical gapless to gapped quantum 
phase transition [7lll0| takes place when /i = 0. 
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FIG. 2: Plots of integrated momentum distribution (in units 
of Likp/n) a) n2,(k^ = 0,fcy) versus fey/fcp b) ny(fcx = 0,fcz) 
versus kz/kp for l/(fc|ap) — —I {+), (hollow squares) and 
1 (solid squares). Insets show: a) n2.{k^ = 0, fcy — 0); b) 
%(fcx — 0,kz— 0) versus l/(fc|ap). 

In Fig. 121 we show the integrated momentum distri- 
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bution ni(k) — J2k '^o(k) = (Li/27r) J dfcino(k) dur- 
ing the evolution from BCS to BEC, where Li is the 
length along i**^ direction. With increasing interaction 
strength, fermion pairs become more tightly bound, and 
thus, ni{\i) becomes broader as fermions with larger mo- 
mentum participate in the formation of bound states. 
Notice also that, while the nodes of the order parame- 
ter are averaged over upon integration, they play an 
important role when either ky or kx is integrated. To be 
specific, at k = 0, riz(O) decreases continuously as a func- 
tion of coupling from BCS to BEC regime, and vanishes 
for l/(fc|ap) -> oo. However, ny{0) = Ly{2M ^lY^^ /n 
vanishes with coupling for /i > in the BCS side, and 
remains zero ny{Q) = for /i < in the BEC side. Thus, 
the qualitative difference between nz(k) and ny(k) [or 
nx(k)] around k = explicitly shows a direct measur- 
able consequence of the gapless to gapped quantum phase 
transition when fi = 0. 

Next we discuss, p-wave superfluidity near Tc. For T = 
Tc (Ao = 0), iVo = X^k "-F[C(k)], corresponds to the num- 
ber of unbound fermions. Here, nF{w) = l/[exp{Pw) + 1] 
is the Fermi distribution. The fluctuation contribution 
-^fluct is obtained as follows. The matrix F~^[q) can be 
simplified to yield 



L- 



l-nF(e+)-nF(C-) 



IVi 



|r(k)p, (6) 



which is the generalization of the s-wave case . Here, 
L-\q) = F-[^{q), and ^± = ^(k ± q/2). The result- 
ing action then leads to the thermodynamic potential 

f^gauss =^0+ ^^fluct, where fJfluct = Y.qMf^I'il)]- 

The branch cut (scattering) contribution flgc to r^fluct 
is obtained by writing [3L{q) in terms of the phase 
shift (5(q, = Arg[/3L(q, + iO+)], leading to fisc — 
-^'^ E^C. nBiw)Sici,w)dw, where u;* = |q|V(4M) - 

2/i and 6{c[,w) = d{q,w) — (5(q, 0). Here, nQ{w) = 
1 / [exp(/3?ii) — 1] is the Bose distribution. For each q, the 
integral only contributes for w > w*, since S{q,w) = 
otherwise. Thus, the branch cut contribution to the num- 
ber equation Ngc — —dflsc/dfJ- is given by 



1 



E 

q 



dnB{w) 



riBiw) 



d_ 
9/i 



S{q,w)dw, (7) 



where w = w + w* . 

When Op < 0, there are no bound states above Tc and 
Nsc represents the correction due to scattering states. On 
the other hand, when Op > 0, there may also be bound 
states in the two-fermion spectrum, represented by poles 
with w < w*. For arbitrary l/(fc|ap), the evaluation of 
the pole (bound state) contribution iVbs requires heavy 
numerics. However in strong coupling. 



(8) 



where Wq = |qp/(4M) and fiB = —Eh + 2/i. Here, we 
use l/i^TTg) = j^u |r(k)|V[2e(k)-£;b] to express Eq. © 
in terms of binding energy < 0. Notice that the ex- 
pression for A'bs given above is good only for couplings 
where /ib < 0. Thus, our results for ko « 200fcF are not 
strictly valid when < l/(fc|,ap) < l/(fc|,a*) ^ 5, where 
a* corresponds to fiB = 0. Therefore, in this region we in- 
terpolate. The binding energy in the BEC regime is i?b = 
—2/{Mkoap) (when k^Up » 1). This result is consistent 
with a T-matrix calculation where = l/(Maprp) 
with rp = -2/(fc2ap) - ^fc2/(4M2l^)Ejr(k)|V62(k). 
This leads to Vp = — fco/2 (when fepOp ^ 1), indicating 
that both approaches produce the same result. 
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FIG. 3: Plots of reduced a) critical temperature Tr — Tc/ep 
and chemical potential fir ~ m/^f (inset), and b) fraction of 
unbound Fq — Nq/N, scattering Fsc = Nac/N, bound i^bs = 
Nhe/N fermions at T = Tc versus l/(fc|ap). 

To obtain the evolution from BCS to BEC, we solve 
numerically the number N — Nq + Nsc + -^bs and or- 
der parameter equations. In Fig. Oa,, we plot T^- = 
Tc/ep and = /x/ep as a function of l/(fcpOp). The 
weak couphng Tc — (8/7r)eF exp[7 — 8/3 -|- 7rfco/(4A:F) — 
7r/(2A;||ap|)] evolves continuously to the dilute Bose gas 
Tc = 27r[2nB/C(3/2)]2/3/MB = O.lSTep in the BEC 
regime, where 7 « 0.577 is the Euler's constant and 
UB = n/2 = A:|/(127r2) is the density and Mb = 2M 
is the mass of the bosons. However, the saddle point 
To « i?b/[21n(£;b/eF)^^^] increases with l/(A:|op), and 
is a measure of the pair dissociation temperature jl3l |. 
Notice that, the ratio of A(A:f)/Tc = Aork(fcF)/Tc in 
the BCS limit is Sir/e'^. The hump in the intermediate 
regime is similar to the one observed in fermion-boson 
model 12]. Furthermore, similar humps were also calcu- 
lated in the s-wave case however, whether they are 
physical or not may require a fully self-consistent numer- 
ical approach. 

The weak coupling = ep evolves continuously to the 
strong couphng /x = —l/{Mkoap) (when fegCp 3> 1) lead- 
ing to fj. — Eh/2. Notice that /i crosses the bottom of the 
band at l/(A:pap) w 0.5, i.e., after the two-body bound 
state threshold l/(fcFap) = is reached. The evolution 
of At at r = (Fig. HJ) and T = Tc (Fig. is similar, 
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but very different from s-wave However, another 

result for fi versus l/(fc|ap) at T = Tc (much hke the 
s-wave case) was obtained in Ref. [l^ using a fermion- 
boson model. In Fig. ^p, we also plot the fractions of 
unbound {Fq — Nq/N), scattering (Fgc = Ngc/N), and 
bound (Fbs = N\js/N) fermions as a function of l/(fcpap). 
While Nq (A^bs) dominates in weak (strong) coupling, A'sc 
is dominant at the intermediate regime. 

Next, we study the evolution of time-dependent 
Ginzburg-Landau (TDGL) equation near Tc. We expand 
the effective action ScS around Aq = to fourth order 
in A{q) [l^, and obtain the TDGL equation 



<i,3> 



2M ' 



A{x) = 0, (9) 



in the real space x = (x, t) representation. The 
time-independent expansion coefficients are given 
by a - 1/(4^.9) - ^|r(k)| V[2C(k)], and 

/32xrfc,fcj/[16AfC(k)]}|r(k)|2, where (5y is 
the Kronecker delta, X — tanh[/3^(k)/2] and 
Y = sech^[/3^(k)/2]. The coefficient of the nonlin- 
ear term is & = {^/[4?^(k)] - l3Y / [8^^ (k)]} \T (k)\^ . 
The time-dependent coefficient has real and imaginary 
parts, and is given by d = Ek ^|r(k)| V[4C^(k)] + 

i/37V(eF)M^/^e(Ai)/(32eoey^): where e(^) is the Heav- 
iside function. As the coupling grows, the coefficient 
of the propagating term (Re[(i]) increases, while the 
damping term (Im[(i]) decreases until it vanishes beyond 
II — 0. This indicates that the dynamics of A(x) is 
undamped for /i < 0. For completeness, we present next 
the asymptotic forms of a, 6, Cij and d. 

In weak coupling {^i — ep), we find a = Ky^ln{T/Tc), 
b = 2K„eFC(3)/(57rT2eo), c^^ = Cyy = = 
7KweFC(3)/(207r2T2), = 0, and d = Kw[l/(4eF) + 
«7r/(8Tc)], where — eF-/V(eF)/(47reo) and ({x) is the 
zcta function. By rescaling the order parameter "^^(x) — 
y^b/ Kv/A{x) one obtains the conventional TDGL equa- 

tion e*w + l^-wP^w - Ez(Cgl)'V2«'w + TGL^t^w = 

with characteristic length £,f^ = Cii/{2Ma) 



and time r 
r)/Tc with 



—id/ a 
« 1, 



tgl/c scale. Here, e 



= (Tc- 
11/3 = 

v/7C(3)/(207r2)(eF/Tc), and tql = -*/(4eF) + vr/(8Te). 
The system is overdamped since Tc ^ ef reflecting 
the presence of two-fermion continuum states into which 
Cooper pairs can decay. 

In strong coupling (eo ^ ^ Tc), we find a — 



9Ks/(2567reo), 



Ks(5y/16, and 



Eb\)/8. b 

where Kg = N {ep) / (iy/ep^) . By rescaling 
the order parameter ^'5(2^) — VdA(x) one obtains the 
conventional Gross-Pitaevskii equation for a dilute gas 
of bosons /ZB*s + C/bI^-sP^s - V2*s/(2A/b) - idt-^s = 
with bosonic chemical potential ^b = —a/d = 2^ — Eh, 
mass Mb = Md/cu — 2M, and repulsive interactions 



Ub = b/d'^ = 187r/(Mfco). In this regime, kpCcL = 
[7rfco/(36/cF)]^^^ is independent of Op, and is infinitely 
large when kg/k-p — > 00. 

The evolution of follows from ("^ql)^ = 

c,,/[2Mrc(aa/9T)] where da/dT = Ek {Y/i-iT^) + 
(9Ai/9T) [y/ (4Tf (k) ) - X/ (2f 2 (k) )] } I F (k) 1 2 . Notice that , 
dfi/dT vanishes in weak coupling, while it plays an im- 
portant role in strong coupling. The evaluation of dfj,/dT 
for intermediate coupling is very difficult, thus an inter- 
polation for connecting the weak and strong coupling 
regimes is shown in Fig. ^p. While representing the 
phase coherence length is large compared to interparticle 
spacing in both BCS and BEG limits, it has a minimum 
in the unitarity region l/(fcpap) 

In the same figure, we 
and the average Gooper 
-(^(k)|V2|^(k))/(^(k)|^(k)), 
A(k)/[2i;(k)] is the T = 



0. 

also 

pair ^pg^j^ 
where V'(k) 
pair 



compare 
size S.I 



zz 
GL 



wave function. 



Notice that ^pair is a decreasing function of interaction 
(while is not). The limiting value of ^pair in strong 
coupling is controlled by kp/ko- Furthermore, ^pair 
has a cusp (non-analiticity) when /i = 0. This cusp is 
associated with the change in i?(k) from gapless (with 
line nodes) in the BCS to fully gapped in the BEG side. 

In conclusion, we analysed the evolution of superfluid 
properties of a 3D dilute p-wave Fermi gas from weak 
(BGS) to strong (BEG) couphng regime as a function of 
scattering volume at temperatures T — and T = Tc. 
We discussed the order parameter, chemical potential, 
average Gooper pair size, and momentum distribution at 
the ground state (T = 0). We also discussed the crit- 
ical temperature Tc, chemical potential and number of 
unbound, scattering and bound fermions at the normal 
state (T — Tc). Lastly, we derived the TDGL equation 
for T w Tc and extracted the GL coherence length. 
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